Anomalous pressure in fluctuating shear flow 
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We investigate how the pressure in fluctuating shear flow depends on the shear rate S and on the 
system size L by studying fluctuating hydrodynamics under shear conditions. We derive anomalous 
forms of the pressure for two limiting values of the dimensionless parameter A = SL 2 /u, where v 
is the kinetic viscosity. In the case A <C 1, the pressure is not an intensive quantity because of the 
influence of the long-range spatial correlations of momentum fluctuations. In the other limit A > 1, 
the long-range correlations are suppressed at large distances, and the pressure is intensive. In this 
case, however, there is the interesting effect that the non-equilibrium correction to the pressure 
is proportional to S 3 ' 2 , which was previously obtained with the projection operator method [K. 
Kawasaki and J. D. Gunton, Phys. Rev. A 8, 2048, (1973)]. 

PACS numbers: 47.27.-i, 05.20.Jj, 05.40.-a 



It has been confirmed over a century of study that 
hydrodynamic equations describe macroscopic flow with 
high accuracy. However, a microscopic foundation of 
fluid mechanics has not yet been established except in 
the dilute limit gas. There seems to be no basis in mi- 
croscopic physical laws on which to establish the validity 
of the local equilibrium assumption inherent in hydrody- 
namic equations. Indeed, it may be the case that any 
proper equation of state must be incorporated with non- 
equilibrium effects. 

The anomalous form of equation of state under non- 
equilibrium conditions was first discussed by Kawasaki 
and Gunton . They derived a non-analytic dependence 
of the pressure tensor on the shear rate for an uniformly 
sheared simple fluid. Subsequently, the same problem 
was studied by several authors in detail Q, 0, HI 0] , and 
the non-analytic response due to mode-coupling effects 
has now become evident. Nevertheless, with regard to the 
normal stress, such predictions have not been fully con- 
firmed by laboratory experiments or by numerical simu- 
lations [5|- 

In addition, quite recently, another type of anomalous 
non-equilibrium pressure has been reported by Aoki and 
Kusnezov 0. They have studied numerically heat con- 
duction problems of anharmonic oscillator models and 
have found that the nonequilibrium correction to the 
pressure is nearly proportional to the system size. This 
finding is remarkable because it implies the absence of the 
intensivity of the pressure in this nonequilibrium system. 
As far as we know, there is no theoretical explanation for 
this result so far. 

In this Letter, we propose a unified understanding of 
the two anomalies for a specific simple case, uniform 
shear flow in a fluid. We establish a general criterion 
which can distinguish between non-intensive and non- 
analytic nature of the pressure. We also give a brief com- 
ment on numerical experiments which disagree with the 
mode-coupling theory. 

The key idea of our study is to establish a possible re- 
lationship between the pressure anomalies and the long- 



range spatial correlation of momentum fluctuations. It 
has been recognized that a lack of detailed balance is re- 
sponsible for various distinctive features. One of them is 
the generic existence of long-range spatial correlations of 
fluctuations of conserved quantities 0. For example, a 
correlation of momentum fluctuations becomes spatially 
long-ranged under shear conditions [j| . 

However, because the correlation function known at 
present is strongly divergent in a long wavelength limit, 
a certain appropriate length scale may be introduced as 
a cut off. Therefore we first need to know the momen- 
tum correlation function over the whole region of length 
scales. It may be suitable to study fluctuating hydrody- 
namics for this purpose. 

Model: Let Vi(r,t) be a fluctuating velocity field in 
an incompressible fluid with a constant temperature T 
that is far from the critical point. The time evolution 
of Vi(r, t) is described by the continuity equation of the 
momentum, 



pd t Vi + djllij = 0, 
where the momentum flux tensor n is given by 



n, 



pviVj +pS i3 - vpidiVj + djVi) + Si 



(1) 



(2) 



Here, v is the kinematic viscosity, p is a constant den- 
sity, and s represents a Gaussian random stress tensor 
satisfying the fluctuation-dissipation relation [T(ij]. 

(s ik (r,t)s lm (r',t')) = 2pTv [(5 t i6 km + SimSki) 



5 3 (r-r')S(t-t').(3) 



Throughout this Letter, the Boltzmann constant is set to 
unity. The auxiliary field p in (J2J) can be replaced by use 
of the incompressible condition divi = 0, which implies 
the relation 

P = Pb - did m Shn - pA^ 1 (div m )(d m vi), (4) 
where p% is a constant. 



2 



We consider the system in a three dimensional space, 
for which — oo < x,z < oo and —L/2 < y < L/2. The 
boundary conditions are imposed satisfying 



v(x, L/2, z) = v(x, —L/2, z) + SLe x 



(5) 



which are chosen so as to make the analysis as simple as 
possible. We study linear fluctuations around the average 
shear flow 



(6) 



(v(x,y,z)) = {Sy, 0,0). 

Without loss of generality, we set S > 0. 

From (JIJ, we have dj (H y j) = for the statistical 
steady state. Because this represents a force balance 
equation, the pressure in the y direction is given by (Tl yy ). 
Then, using J2J and Q , we write the pressure p y as 



P y = PB +p{vy) - pA 1 {{div m )(d m vi)) 



(7) 



Now, the problem is to determine how p y depends on the 
shear rate S and on the system size L. 

The most general form of fluctuating hydrodynamics 
consists of a set of stochastic evolution equations for the 
energy, density and momenta |l0j |. In our model, the 
thermal diffusion constant and the sound damping con- 
stant are assumed to be much larger than the kinematic 
viscosity. This assumption is made for simplicity here, 
otherwise the calculation of the stress tensor becomes 
more complicated 0, • 

Now, we make a comment on the local equilibrium as- 
sumption that is involved in hydrodynamic descriptions. 
Let £ be a microscopic length, such as a mean free path, 
and define a dimensionless parameter e as 



es 



(8) 



In order to investigate the physical meaning of this pa- 
rameter, we first consider the case e — > 0. In this case, 
the local equilibrium assumption is valid because of the 
separation of scales. This implies that the shear flow 
does not influence the thermodynamic properties of the 
system. However, if e becomes of order unity, the molec- 
ular motion is violently disturbed by the shear. Although 
it may be interesting to study such systems, we cannot 
use fluctuating hydrodynamics for that aim. From this 
brief consideration of these two limiting cases, we regard 
the parameter e as representing the extent of departure 
from the condition of local equilibrium states. Because 
we seek the non-equilibrium correction to the pressure 
under shear within the framework of fluctuating hydro- 
dynamics, we consider the case of a small but finite e. 

As the final condition on our model, we assume L 3> i, 
so that the continuum description is valid. Note that our 
model does not include the parameter I explicitly. It is 
included only as a cut-off length when we encounter an 
ultraviolet divergence. 



Dimensional analysis : There are two independent di- 
mensionless parameters in the system: 



, L 2 S T 
A = , D 



v pv 2 L 
Quite generally, we can write 

T . / L 2 S T 



Py~Py q =nf 



L 3 \ v pv 2 L 



(9) 



(10) 



where p y q is the equilibrium pressure. We have assumed 
that p y — p y q does not have ultraviolet divergences. The 
validity of this assumption is not obvious, and it should 
be confirmed by a concrete calculation. 

Because we carry out a linear analysis of the fluctua- 
tions, we replace f(X,D) by / (A) = f(\,D = 0). Then 
the problem is now reduced to deriving the form of the 
function /o(A). However, because it is still difficult to 
find a general solution even in this simplified problem, 
we study two specific cases, A< 1 and A > 1. 

When A <C 1, all quantities may be expanded in powers 
of A. Due to the reflection symmetry, we have /o(A) ~ 
Co A 2 . Here cq is a constant whose value is calculated to 
be 1/11527T below. We thus obtain 



Py - Py q = c a TL 



S 



(11) 



Note that this L dependence, which implies the break- 
down of the intensive nature, is compatible with the re- 
sult by Aoki and Kusnezov |6j] . At the end of this Let- 
ter, we discuss that this anomalous behavior can be at- 
tributed to the long-range correlation of momentum fluc- 
tuations. 

Next, we shall consider the opposite case, A ^> 1. 
Note that this condition does not necessarily imply the 
strong shear. For example, this asymptotic case is re- 
alized by taking the limit L — > oo, while keeping the 
shear rate S small, in which case the shear stress obvi- 
ously remains small. Note also that planar Couette flow 
is linearly stable for all values of A. If we encounter no in- 
frared divergence in the calculation of p y —p y q , we obtain 
/o(A) = ci\ z / 2 , because p y — p y q should be independent 
of L. In this case we have 



Py -P 



cq 



ciT 



3/2 



(12) 



This S dependence is the same as that obtained by 
Kawasaki and Gunton. As we demonstrate below in- 
deed, no infrared divergence is accompanied in the cal- 
culation of p y —p y q , and that calculation yields l|12f> with 
d = 1.06 x 10~ 2 . 

Technical details : We first consider the case A — > oo 
and begin by performing the Fourier expansion 



(r) 



dk 



Ahr 



Vy(k). 



(13) 
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In a statistically steady state, an equal-time correlation 
function can be defined as 



(Vy(k)Vy(k')) = Cyy(k)S 3 (k + k') 

Equation J7J) is then rewritten in the form 

dk 



p y =PB+ p 



-Cyy(k), 



(14) 



(15) 



where we have used the result ((div m )(d m vi)} = 0, which 
is readily obtained by divi = and ©. From (JTJ and 
@, we derive the equation of the correlation function 
Cyy(k): 



dCyy(k) 

dk v 



2 vk 2 



h h 



Cyy(k) 



2vT 



(k x + k z ). 



(16) 



We can solve this equation by introducing the new 
wavevector k' — k + Stk x e y , where t is a fictitious time 
parameter. Using the standard manipulation |H, lllj. we 
obtain 



Cyy(k) — 



2vT 



P 



(k 2 + k 2 ) / ( ft e -2vk 2 (t+t 2 Sk x k y +t 3 S 2 kl/3) 



(l + 2tSk x k y + t z S z k 2 y. (17) 

Combining (|15l) and (|17f) . we find an integral expres- 
sion of the pressure [l^- When S = 0, this merely gives 



the equilibrium pressure p y q , which can be written as 



P? 



PB 



bTt 



(18) 



Here I is introduced as a short scale cut-off to avoid the 
ultraviolet divergence, and & is a numerical constant. Be- 
cause pb remains undetermined, this ultraviolet diver- 
gence is renormalized into it appropriately so that p y q 
coincides with the observable equilibrium pressure. In a 
nonequilibrium case (S > 0), we find that the correction 
to the equilibrium pressure is actually expressed in the 
form (13), and C\ is given by 



ci 



1~tF l f°° l 



with 



g(t) 



dOA(O)- 



B(Q) + tC(n) 



,3/2 • 



(19) 



(20) 



[l + £B(Q) + f C(O)]' 

Here Q = {0,4>), dfl = d<j)d6 sin6 and A, B, C are 
polynomials of sin#, cos0 and sin0, which are deter- 
mined from (|17|l . From numerical integrations, we obtain 
ci = 1.06 x 10^ 2 . 

In the case A < 1, we can utilize the calculation per- 
formed above by simply replacing J dk y (2Tr)~ 1 e lk y y ■ ■ ■ 
with J2 n L~ 1 e l27 ' ny / L ■ ■ ■ in l(T3|l . In this way, an expres- 
sion similar to (|17f) with the replacement k y — > 2n/L is 
found. Then the expansion of this correlation function in 
powers of A leads to the desired form of the nonequilib- 
rium pressure in (|llfl . whose numerical factor is obtained 
as c = 1/11527TW 2.8 x 10" 4 . 



Long-range correlations : Here we demonstrate that 
the anomalous forms of the pressure (|TT1) and (|T2|) are 
closely related to the long-range correlations of momen- 
tum fluctuations. For simplicity, we restrict our attention 
to the fluctuations with k = (fc,0,0). Two asymptotic 
forms of the correlation function are given as, 



c„„(fc,o,o) 



T 
P 



1 S 2 



2 v 2 k 4 



for fc 2 > S/u, and 



C yy (k, 0,0) 



T 
P 



1/3 



u 2/3 k i/3 ' 



(21) 



(22) 



for k 2 <C S/u. Comparing these expressions with 111|) 
and l|12|l . the dependence of the pressure on the long- 
range fluctuations in each case becomes clear. The 
1/fc 4 dependence in l|21|l reflects highly anomalous be- 
havior of the fluctuations with small wavenumbers |13| . 
Note that there is a range of wavenumbers that satisfy 
k 2 3> S/v but are still much smaller than the charac- 
teristic wavenumber of equilibrium density fluctuations. 
On the other hand, \\1'2\\ shows that this long-range cor- 
relation is suppressed at scales larger than I = yjvjs, 
crossing over to a weaker correlation (See also Fig. QJ. 
This is equivalent to the stronger power law decay of 
the correlation function than 1/r in a real space [lj| . A 
new length scale I characterizes this crossover, which is 
intrinsic in the nonequilibrium system considering now. 



F(x) 10 




FIG. 1: Plots of numerically evaluated scaling function F(x), 
where the correlation function C yy (k) is written in the form 
C vy (k) = p' 1 T(l + F{kl)). Note that p~ 1 TF{kl) corresponds 
to the nonequilibrium correction to the momentum correla- 
tion function. The solid line and the dotted line represent the 
asymptotic functions calculated from the exact integral form 
of F(x) for large x S> 1 and for small i < 1, respectively. 
These asymptotic functions are equivalent to the results ll'H 
and 1221 . 

These results provide the following physical picture. 
The momentum fluctuations exhibit the long-range cor- 
relation described by J2U and (|2"2")l . When L <C I, this 
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correlation yields the non-intensive contribution to the 
pressure given in Ijlll) . On the other hand, when L is 
chosen to be sufficiently large, the long-range correlation 
is suppressed at scales larger than I, and this leads to the 
non-analytic shear rate dependence given by l|12|l instead. 

Discussion: In a nonequilibrium system, an external 
field having a spatial gradient (i.e. shear flow) induces 
the coupling of fluctuations with different wavevectors. 
In a perturbative expansion to the lowest order in shear 
rate, a correlation function has the same form as (|21[) . 
This may cause an infrared divergence in the calculation 
of p y —p y q - Roughly speaking, (f[Tj) is obtained when the 
cut-off scale is chosen as L, while (|12ll is obtained when it 
is chosen as I. We have demonstrated the validity of this 
intuitive argument by computing the correlation function 
rigorously. 

The predicted finite-size dependence of the pressure 
under the condition A <C 1 is fairly striking. Since 
/ = \J v I S may be less than 1 mm for water at standard 
temperature and pressure and for experimentally acces- 
sible shear rate, it is possible to design an experimental 
device corresponding to the the condition A < 1. On 
the other hand, from ©, dJ and JTBJ), (p y - p^) / p c y q is 
found to be proportional to e 3 / 2 in the case A > 1, and e 2 
in the case A <C 1, respectively. This indicates that the 
nonequilibrium correction vanishes in local equilibrium 
states (i.e., e — > 0), as expected. For example, e is less 
than 10 -8 when the shear rate is 10 3 sec -1 for water con- 
sidered above. Thus, unfortunately, the nonequilibrium 
correction may be too small to observe for simple fluids 
under ordinary experimental conditions in either case. 

This fact, however, does not eliminate the fundamen- 
tal significance of the present study. In constructing a 
theoretical framework of non-equilibrium statistical me- 
chanics, this 'e-effect' must be taken into consideration, 
because a statistical distribution can be reduced to a lo- 



cal canonical ensemble only in the limiting case e — > 0. 
Furthermore, making use of the experimental technique 
developed in the micro- and nano-fluid studies may en- 
able us to detect the novel nonequilibrium effect for real 
fluids in the near future. We also expect that our findings 
stimulates further numerical studies such as molecular 
dynamic simulations. 

Recently, Marcelli et al. reported the analytic depen- 
dence of pressure in shear flow observed in simulations 
of non-equilibrium molecular dynamics Because the 
shear conditions were chosen to satisfy e > 1 and A 3> 1 in 
their simulations, their result is not directly comparable 
to ours. However, we believe that the qualitative feature 
of our result does not change even in the case e > 1. We 
suspect that their model system does not exhibit long- 
range momentum correlation because of the absence of 
the local conservation of momentum. That might be the 
reason why their result is inconsistent with the existing 
mode-coupling theories P, 0, Q) and ours. We expect 
that measuring momentum correlation functions will help 
to resolve the discrepancy, probably favorably with the 
mode-coupling theory. 

Finally, we give a general remark on the possibility 
of describing an equation of state for a non-equilibrium 
steady state using a thermodynamic function. An ex- 
perimental test of this possibility has recently been pro- 
posed in |15| . in which the intensivity of the pressure is 
postulated. Therefore, the recovery of the intensivity in 
our calculation has a significant meaning with regard to 
the construction of thermodynamics extended to steady 
states. 
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